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We prove that QCD in the e-regime of chiral Perturbation Theory is equivalent to chiral Random Matrix 

, Theory for zero and both non-zero real and imaginary chemical potential /i. To this aim we prove a theorem 

^ ■ that relates integrals over fermionic and bosonic variables to super-Hermitian or super-Unitary groups also 

\jQ ' called superbosonisation. Our findings extend previous results for the equivalence of the partition functions, 

, spectral densities and the quenched two-point densities. We can show that all fc-point density correlation 

CO ' functions agree in both theories for an arbitrary number of quark flavours, for either jj, — or jj, taking 

' real or imaginary values. This implies the equivalence for all individual fc-th eigenvalue distributions which 

. are particularly useful to determine low energy constants from Lattice QCD with chiral fermions. 
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1 Introduction 



The application of chiral Random Matrix Theory (ehRMT) to QCD, initiaUy introduced at zero chemical 
potential in [Tj, has led to many exact analytical results for the Dirac operator spectrum. These have in turn 
found applications in comparison to Lattice Gauge theory by various groups (for reviews see The deeper 

reason for this analytical understanding of the Dirac spectrum in the phase with broken chiral symmetry can be 
traced back to the relation between chRMT and the epsilon regime of chiral Perturbation Theory (exPT). Being 
based on global symmetries alone there are three classes of chiral symmetry breaking patterns, corresponding 
to the three possible anti- unitary symmetries that lead to the unitary, orthogonal and symplectic chRMT [i]. 
Here, we will only report on the first class containing QCD, which is technically the simplest. Relatively little 
is known about the equivalence to exPT for gauge theories in the two other classes, see however ^5ii6j. 

In this paper we will show that any correlation function of Dirac operator eigenvalues computed starting 
from exPT is the same as the corresponding quantity computed from chRMT in the large- A'^ limit of infinite 
matrices, corresponding to the infinite volume limit of xPT. This equivalence holds both with and without 
chemical potential, being of real or imaginary type. The chRMT we use is the one introduced by Shuryak and 
Verbaarschot [T] for ^ = 0, by Osborn for real n ^ 0, and by Akemann et al. [5] for imaginary /i where the 
two latter are two-matrix models. The one-matrix model for real ^ introduced earlier by Stephanov |9j can 
also be mapped to exPT and is thus equivalent to [Tj . Being technically more involved we don't give details of 
our analogous equivalence proof here. 

The exPT - chRMT equivalence is established by relating the generating functionals of resolvents of an 
arbitrary finite degree, with any number of flavours Nf. In exPT these generating functionals are obtained by 
adding additional fermionic and bosonic quarks, leading to supergroup integrals. For real 7^ it is nontrivial 
to make them convergent, and our result for the corresponding exPT is new. These group integrals are then 
shown to be equal to the corresponding quantities in chRMT at large-A^, given by the expectation value of 
ratios of Dirac operator determinants (also called characteristic polynomials in the chRMT language) . 

The technical tool we use is a theorem we prove also called superbosonisation, which was obtained indepen- 
dently by [ini E] . After writing the ratios of determinants in terms of superfermionic variables and integrating 
out the Gaussian random matrices the theorem allows us to map these expressions to the supergroup integrals 
of e^PT. We only take the large- limit at the very end, without performing any saddle point approximations. 
Let us emphasize that although there are many applications of RMT in Physics [T^] , only in few cases an exact 
map to the underlying microscopic theory has been achieved. 

In chRMT there are alternative ways to compute arbitrary fc-point density correlations, or individual eigen- 
value correlation functions, some without introducing resolvents. One example is the method of orthogonal 
polynomials T3]. On the other hand the supergroup integrals in e^PT are increasingly hard to calculate for 
k > 1. By establishing an equivalence for all generating functionals we can therefore match with all known 
chRMT correlation functions, were they obtained using resolvents or not. For fi — all generating functionals 
in chRMT were computed in [TH [H] . For real /i ^ all building blocks for characteristic polynomials in chRMT 
in the large- limit follow from [16j and |17l . At imaginary /i ^ so far only the bi-orthogonal polynomial 
method has been used in chRMT [8]. 

The explicit results for the fc-point density correlation functions can be most conveniently read off directly 
from [TH] for massless, and from [TH |^ for massive correlations at fj, = 0, for both massive and massless 
correlations from [THI at real /i 7^ 0, and likewise from [5] for imaginary /i. The corresponding individual 
eigenvalue correlation functions that require the knowledge of all density correlations were obtained in 21j for 
/X = 0, in |22j for real fj, and in ^23j for imaginary fi. 

The following results were previously known about the exPT - chRMT equivalence: partition functions with 
only fermions were shown to agree for zero [T] and non-zero fi [241 1251 1261 [3] . The equivalence for partition 
functions with only bosons at /x = follows from ^27j. The generating functional of a {Nf -I- 1|1) supergroup 
integral leading to the spectral density was computed in [2SJ [21] including Nf massless fermions. This work 
was extended to include the quenched two-point density in |30j . in this case too only for /z = 0. Furthermore, 
the quenched density at real /i ^ was computed from both e^PT [14^ and chRMT [7] and found to be in 
agreement. For imaginary isospin /i 7^ the equivalence was established up to the two-point function in [311 15]. 
It was pointed out 32J how in principle to compute the distributions of the fc-th individual eigenvalue from 
exPT, using all fc-point density correlation functions, in order to reproduce previous chRMT results [2T\ . The 
same strategy can be applied for non-zero n [55] [32] . 
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Our equivalence proof for all fc-point density correlation functions thus, in particular, fills this gap. Indi- 
vidual eigenvalue distributions have been used in several groups to compare with Lattice results with exact or 
approximate chiral symmetry, leading to a well defined gauge field topology. 

A further remark concerns universality. On the chRMT side it is well known for /z = that in the microscopic 
large- TV limit the results hold for a much wider class than Gaussian weight functions [HI [33l US] ■ For that reason 
we can restrict ourselves to establish the equivalence for Gaussian chRMT. For fJ, less is known. Both the 
model by Stephanov [5] and by Osborn [7] are Gaussian matrix models. As already mentioned above the 
fermionic partition functions [211 [13 [HI [3] and the quenched density [TH [7] were shown to agree among both 
matrix models, and with e^PT. We can now extend this equivalence to all correlation functions. This agreement 
suggests universality for chRMT in the broken phase to hold at /i ^ as well. 

Finally let us mention that for non-chiral RMT the equivalence to an effective theory for spontaneous flavour 
symmetry breaking in 3D is simpler and has already been fully established at /i = by Szabo [34 . We expect 
that an extension to /i 7^ can be done following the same lines as here. 

Our paper is organised as follows: In section [D the resolvent method is briefly recalled both for Hermitian 
operators (/Lt = 0, pure imaginary) and non-Hermitian ones (/i real), in order to clarify how the A:-point 
correlation functions are obtained from a partially quenched (pq) theory. In sections [3l [4| and [5| the proofs 
of the equivalence for the pq theories for /i = 0, for imaginary fj. and real /i are given, respectively. Here we 
write the averages over ratios of Dirac determinants in chRMT as supervectors integrals, perform the chRMT 
average and apply our theorem relating to exPT supergroup integrals. In particular the form of pq e^PT with 
real /.t ^ is given. The superbosonisation theorem is then presented and proved in section [Hj More technical 
remarks on delta functions and on some integrals used are deferred to appendices [3 and [HI 

2 The resolvent method 

The resolvent method is a way to compute the fc-point eigenvalue correlation function of an operator distributed 
according to a given ensemble (in our case these ensembles are QCD or chRMT). Starting from the expectation 
value of a ratio of characteristic polynomials, this quantity may be considered as the partition function of a 
theory with an additional number of fermionic and bosonic quarks, called partially quenched QCD (or chRMT). 
Both supersymmetric and replica methods rely on the resolvent method. 

In this section we briefly summarise it's idea for the reader's convenience, referring to the literature for 
more detailed explanations [35]. Let us emphasize that the resolvent method is substantially different whether 
one considers a theory with an Hermitian or non-Hermitian operator. QCD with zero or imaginary chemical 
potential is a theory with an anti-Hermitian Dirac operator, the correspondence with an Hermitian operator is 
trivial. QCD with real chemical potential is a non-Hermitian theory. 

We start by showing the resolvent method for Hermitian theories. The fc-point correlation function is 
defined as the expectation value of the product of fc J-functions: 



where the expectation value is computed according to the ensemble of the operator D whose eigenvalues are 
A. Alternative definitions of this quantity may be given, in order to avoid the contact terms obtained through 
delta functions at coinciding arguments r]i — rjj [12| . 

For simplicity wc consider first the easiest case, the spectral density (fc = 1). We define its resolvent as 



Here the argument z is taken to lie outside the support oi D <Z TZ. Given that the density is z-independent 
this integral equation can be inverted as follows |36j : 




(1) 




(2) 



(3) 
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taking the discontinuity of the resolvent when crossing the support. An essential ingredient for this result is the 
fact that all the eigenvalues of an (anti-)Hermitian operator lie on the real (imaginary) axis. This requirement 
is not fulfilled by theories with a non-Hermitian Dirac operator, invalidating these equations. The next step is 
to generate the resolvent ([5]) as a ratio of determinants: 



Gi{z) 



Tr 



z- D 



d /Vet[D-z'] 



dz' \ Vet [D - z] 



(4) 



The generalisation to the /c-point correlation function is conceptually not different, with a bit more attention 
needed in subtracting eventual higher order poles in z in ^11^=1 X^Aee v F^a)- proved that the k- 

point correlation function may be generated through the theory with k additional bosonic and fermionic quarks. 
The corresponding resolvents are defined as 



Gfe(zi, 



,Zk) 



n E 

\j=l AGC.V. 



1 



and can be generated as in eq. ^ 
Gk{zi, 




1 



1 



-/5fc(Al 




J-i Vet\D- Zi] 

j=i J' 



,Afc) . 



(5) 



(6) 



The density correlations are then given by the discontinuities with respect to all arguments (see e.g. 

1 



Pk{m^ ■■■,11k) 



lim 



E n '^J ^'^ ('^i ~ ■ ■■ , Vk- ic^ks) ■ 



(7) 



{cr},(Tj=±l V = l 



In non-Hermitian theories eigenvalues are no more localised along the real (imaginary) axis, and, as a 
consequence, the correlation functions cannot be computed through the discontinuity of the resolvent along the 
support of the eigenvalues. The equation below will be used to generate the two-dimensional (5-function in the 
complex plane instead: 



1 



= 7r(5^(z) . 



(8) 



Because the average in eq. ([2|) is now over complex eigenvalues the resolvent is well-defined everywhere in the 
complex plane, both in- and outside the support of D. From eq. ^ it follows that for the density we have 



pi{z) = -d:,>Gi{z) 

TT 



(9) 



In other words the resolvent is holomorphic outside the support, and non-holomorphic inside. As before the 
resolvent may be generated by differentiating the expectation value of a ratio of determinants as given in eq. 
([1]). However, this case has an additional problem due to the fact that the microscopic field theory (such as 
exPT) requires the introduction of additional bosonic and fermionic conjugated quarks. This phenomenon is 
due to the Hermiticity requirement [371 1241 117] for bosonic field theories to be written in terms of convergent 
integrals, and we will discuss that in more detail in sect. 15.11 

As before the above can be generalised to any fc-point function, leading to the introduction of k additional 
bosonic and fermionic quarks (or to k couples of conjugated bosonic and fermionic quarks, if one needs to 
consider the microscopic theory), with 



Pk{zi, ...,zk)= J]^ Gk{zi, 



■,Zk) ■ 



(10) 
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Let us mention an alternative method to generate all density correlations - apart from non-resolvent related 
ones: it generates the two-dimensional 5-function using the identity 



Sliz) = lim 



K^o 7r(z • z* + k?Y 



(11) 



and hence when considering the sum of the delta functions over all the eigenvalues A of a non-Hcrmitian matrix 
D we have (see [55]): 

<5^(z-A) = lim irr^ 



Aee.v. 



{{z - D){z - D)"^ + k'^Y 



(12) 



3 Zero chemical potential 

We start from the simplest case: the Hermitian chRMT corresponding to QCD with zero chemical potential. 
In principle this theory could be seen as a particular case of the imaginary jj, case studied in sect. IH setting 
= 0. Since computations are easier here, we present this case separately for pedagogic reasons. 
The idea of our approach to prove the equivalence of all spectral properties in the two theories is the 
following: in both theories the resolvent method allows in principle to compute all fc-point correlation function. 
So without taking the discontinuities, or anti-holomorphic derivatives, it is sufficient to prove that, up to an 
irrelevant constant, the integrals describing these generating functions of the resolvents coincide. From this 
knowledge follows that all the k-point correlation functions are equal, and, hence all spectral properties agree. 
In order to actually compute any given correlation function we may thus choose either theory, and within that 
theory we may even use any other equivalent method that is simpler. 



Hermitian chRMT is a very efficient way to describe the spectral properties of QCD without chemical 
potential [HI!]. The equivalent of the QCD Dirac operator in chRMT is given by 



lA 



/-I, 



(13) 



where A is a complex iV+ x iV_ random matrix. It describes a fixed topological sector of QCD with a topological 
charge v = iV_|- — where v will be kept fixed while considering the A''+ -I- N- = N ^ oo limit. 

Thanks to the universality of the microscopic limit of chRMT [181 [33 HI] we can choose a Gaussian weight 
function for the random matrix A. Since we are interested in finding the spectrum and the fc-point correlation 
function using the resolvent method [39l [28j [29l [34] , according to what was said in sect. [2|we consider the more 
general theory with n/ fermionic quarks and rib bosonic ones and compute it's partition function: 



_ Vet[Vf+mf] ' 



= / dAi 



jNTr\A'< A\ 



n7 Vet 


mflN+ 

iA"! 


iA 
mflN__ 




iA^ 


iA 



(14) 



Here dA is the flat measure in the independent entries of the complex matrix A. The quantity defined above is 
usually called partially-quenched (pq) partition function. 

In the following it will be required that Re{mb) > OVfe in order to have convergent integrals generating the 
resolvents in the upper half-plane. The result for Re{mb) < may be easily recovered from the equation abov^. 

At this point we follow the common procedure [12j of writing the determinants in the numerator in terms 
of Gaussian fermionic integrals, and the ones in the denominator in terms of bosonic ones. We introduce two 
sets of and complex-supervector^ in {nb\nf). In this paper we will use the boson- fermion convention 



^In the literature often a diagonal matrix S with elements = ±1 is introduced to be able to work in the whole complex plane 
at once. For simplicity we omitted such a notation. 

^The anticommuting number part of supervectors and its conjugate are independent Grassmann variables. 
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for ordering elements in supervectors. We write these vectors in a matrix form ipg^a, (t'g,i3j "^g a ^^'^ 't'*g /3 where 
Latin indices run over the (nb|n/) superflavours and Greek indices run over the (N-) eigenvalues, 



E 

g=-nf 



(15) 



Here V'*; is a shorthand notation for the product of the flat measures of the independent entries of 

the supervectors. The integral above depends on the random matrices only in a Gaussian way: 



Exp [~aN Ao,,pAl p - lA^^p {rgo.<l^g,p) - ^A^ p {4'l0iJgo)] ■ 
We can thus perform the Gaussian integration by completing the squares: 



(16) 



Zpq cx / (i('(/',-0*»*) Exp 



-Str 



g,/3 



X Exp 



-Str 



(17) 



In the following we introduce the matrix notation Mgh = N ■ Sghmg, anticipating the correct scaling later in 
the large- A'^ limit to obtain A^-independent quantities in the thermodynamic limit [41 . This expression is a 
Gaussian integral in any of the two sets of supervectors. We can easily perform one supervector Gaussian 
integration (0, (p^) obtaining a superdeterminant as a result: 



Zp, cx / d(V',V'*)e-^*'^K-2:„^9.-®^^,J Sdet 



-JV_ 



(18) 



The equation above depends on the supervectors only trough the sum of external products. We can now 
apply the superbosonisation theorem in eq. (j55p . 



1 

^ ^9,c. «) V'I,a ^ Ugh e Gl{nb\nf) 



(19) 



where Gl{nb\nf) is the maximal Riemannian submanifold of the linear group in the {nh\nf) superspace |42j : 
matrix representations of this manifold are possible, their elements are supermatrices with Gl{ni,) /U (nf,)©C/(n/) 
as a base manifolcjfl and 2nb ■ nf Grassmann variables. The Haar measure on this manifold is called dfj-nici) 
and may be expressed in terms of usual integrations trough (see eq. (IMl) ): 



dfiH 



Gl(nb\nf) 



dHid{Hi) / dixu{H2)Vet[H2T' j dOdQ^Sdet 



Hi ef 
e H2 



nf-nt, 



(20) 



Here, dHi is the flat measure on the Hermitian matrices, 9{Hi) is the product of the step function in the 
eigenvalues {9{Hi) > <;=4» Hi is positive definite), d/ijj is the Haar measure on unitary matrices H2 and 
dO, dO^ is the flat Grassmannian integration in the independent entries of the boson-fermion blocllf|. Other 
parametrisations of Gl{rn,\nf) have been provided for some specific values of (nf,|n/) [521 1331 EH] ■ 
The result is {N+ = v + N^)\ 



cx / d^jiH{U) Sdet[U]"+ e 



N+ ^-Str[^Mu] 



Sdet 



n -AT. 



u 



Gl{nb\nf) 



dn„{U) Sdet [[/] " e"^*" *^^] Sdet 



1 + a—M-U' 



(21) 



■^The bodies, that are the complex number parts of the commuting numbers, of the boson-boson block and fermion-fermion 
block belong to a manifold called base manifold. 

^The Grassmann variables j and @t are independent real Grassmann variables. 
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All the equations above hold for finite N. As a last step we can perform the N ^ oo limit. Here we send 
the masses to^ to zero while keeping Mgh fixed. Therefore the N oo merely leads to an expansion of the 
superdeterminant, without the need to perform any other approximation: 



lim Zpg = [ diiH{U)Sdet[UY e" 

'-^"^ JGUn^Uit) 



V Str\-\MU-aM-U' 

d^lHiU)Sdet[U]''er^^''^i^'■''+"■''"^ . (22) 



OC 

lGl(nb\nf) 

This equation is equivalent to the partially-quenched partition function e.g. in [29], after matching parameters 

miNJ-^nii—. (23) 

We denote by V the volume in xPT and by S the chiral condensate, the first low energy constant in xPT. 
Applying the resolvent method to both eq. in the limit N — > oo, and eq. (|22p while keeping eq. finite, 
we obtain the claimed e^PT-chRMT equivalence at /i = of all fc-point correlation functions in the microscopic 
limit. 



4 Imaginary chemical potential 

The simplest generalisation of the arguments above is obtained by adding imaginary chemical potentials. This 
is a Hermitian theory too and, hence, the resolvent method may be applied in the very same way as before. The 
only difference with the former section lies in a different random matrix Dirac operator, containing one more 
random matrix. This model introduced in [8] follows the idea of [7] in assuming that the chemical potential 
term in non-diagonal in matrix space. This apparent complication by adding more random variables in fact 
makes the model simpler: for two different chemical potentials one can go to an eigenvalue basis and use 
bi-orthogonal polynomials. Below we don't need to diagonalise the matrices as we explicitly integrate them 
out. The equivalence between chRMT and exPT we show thus holds for any number of different chemical 
potentials. However, the integrals and density correlations are so far known explicitly only when having two 
different chemical potentials for any number of flavours. 

In principle one could also follows the idea of Stephanov [S] with only one random matrix and the /i-term 
proportional to unity. However, for technical reasons the computation is more involved and although we don't 
show it here it leads to the same result. 

The Dirac operator for a quark with mass m f and imaginary chemical potential is given by [8] : 



Vf + JTlf 



TTif 1 



+ iHf 



St 



lA + i^fB 

m/lAr„ 



(24) 



and likewise for boson masses nib and chemical potential /if,. Here A and B are complex x 7V_ random 
matrices with Gaussian weights. The issue of universality is more subtle here because the matrices A and B 
will couple after changing variables. We refer to [5] for a more detailed discussion. 

We will prove the equivalence of the partition functions for a generic theory with rih bosons and nj fermions. 
The equivalence of the spectra will then follow, applying the resolvent method. As before we will consider only 
bosonic masses with positive real parts, 



Zpq — 



'YT/VetiVf + mf]' 
nr T^et [Vh + nib] 



Vet 


niflN+ 
iA'' + ififB'' 


lA + ififB 
m/l7v_ 


WVet 


niblN+ 
iA^ + iufB^i 


iA + ifibB 
niblN- 



(25) 



From this point on, most of the steps are equivalent to the ones performed in the previous section: we 
introduce two sets of iV_|_ and complex-supervectors in {nb\nf) to write the ratio of determinants as a 
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Gaussian integral 



J dAdB d(V',V'*,0,0*)Exp [-aN {Tr [A^ A] +Tr [B^ B])] 



(26) 



X Exp 



E 

9=-"/ 



V', 



9,o> 
9,/3 



We reduce all the dependence on random matrices in terms of Gaussian functions: 

Exp [^aN A^^pAl p ~ lA^j, (V'^a'/'s,/?) " (<^s,/3^s")] 

X Exp {-aN Ba^fiB*^j3 - ifigB^^o (V'ga^g,/?) - if^gK^p (^g./sV'sa)] ■ 

Next we can perform the Gaussian integration completing the squares: 



(27) 



Zpg (X d{'4>,'4)*,4),4)*) Exp 



-Str 



X Exp 



X Exp 



-Str 



aN 



^Str 

aN 



t 



(28) 



Here we introduce again the matrix notation Mgh = N -SghTrig, as well as ^^^^ = ^J^^gh^J^g^ anticipating also the 
proper scaling of /i with N below. This expression is a Gaussian integral in any of the two sets of supervectors. 
As before we explicitly integrate one set of supervectors and express the remaining external product in terms 
of an integration over Gl{n})\nf). The result is: 

-JV_ 



N aN 

We can use again the theorem of section [5] obtaining 



Zpq CX 



diXH{U)Sdet[UY e ^'''\.~^''^\ Sdet 



1 + cr-^Af • U-^ + ^B^f'^UB^f'^U-^ 

I\ -i- 1\ 



-N- 



(29) 



(30) 



' Gl{nb\nf) 

This result is again exact for any finite N. If we now take the large- limit while keeping M and S'^-' fixecH 
we obtain finally: 



N 



lim Zpq = 



' Gl(nt\nf) 

where our only approximation has been the expansion of the superdeterminant. 

This expression is equivalent to the e^PT effective partition [31j , where in order to match we use the following 
relations [^[5] 



V 2 2 ' 

M.,Va ^^1D,F^^ . 



(31) 



Here /^d is a dimensional constant, instead of the ^ used above that is dimensionless. The additional parameter 
-FV is the pion decay constant, the second low energy constant in the leading order chiral Lagrangian. Thus 
we have established the exPT-chRMT equivalence for all fc-point correlation functions with imaginary chemical 
potential. 



^In the theory with real chemical potential in the next section this limit is called weak non-Hermiticity limit [44) . While we 
inherit the same scaling here our operators are always Hermitian. 
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5 Real chemical potential 



Before starting the proof for real chemical potential it's better to spend a few words on the definition of theories 
with bosonic quarks. It may help understanding both the result and the proof. 

5.1 Bosonic partition function with /tt 7^ 

The sign problem makes it impossible to write the inverse Dirac operator as a convergent bosonic integral: 
Hermitisation is needed [151 and a new particle content of the theory comes out [371 IHl HZ] • 

Let's consider the simplest example, = 1, nf — 0, following 17J. The partition function is given by: 

,Vet [D + 11^0+171]) ' 

In the microscopic theory (QCD) the inverse power of the determinant should be drawn from the integration of 
bosonic quarks. This integration is well defined only if the integral is convergent, that is, if the Hermitian part 
of the operator is positive (or negative) defined; this requirement is fulfilled when Dirac operator is given by a 
non-Hermitian (kinetic) part plus a mass term (proportional to the identity operator), but this is not the case 
for a Dirac operator suffering from the sign problem. 

The way to bypass this problem is to add an additional bosonic conjugate quark and a fermionic one with 
mass — m* and chemical potential — /i, respectively, 

1 \ _ / Vet[D- fija-m*] 



Vet [D + /170 + m] / \ Vet [D + ^70 + m] Vet [D - ^70 - 

/ Vet [D ~ fijQ — TO*] 



I Pet 



D + /170 

D — /i7o — TO* 



(33) 



(34) 



The operator in the denominator is now anti-Hermitian and hence has only pure imaginary eigenvalues. If we 
consider the matrix 

e D + fijo + 

D — /i7o — TO* e 

(sending e ^ at the end) we have an operator fulfilling the positivity requirements to have a properly defined 
theory. This process is usually called Hermitisation ^45> i46j. and although it could seem as just a mathematical 
trick this is not the case: the possibility of writing the inverse determinants as convergent integrals is a feature 
of the underlying theory. This procedure modifies the form of the mass term in e^PT [24[ I17j , and for the most 
general result we refer to eq. (|49|) , where the generalisation to any number ni, of bosonic uncoupled quarks is 
straightforward: (iVf,|0) {Nb + Nj^\N^). Despite this regularisation is a feature of the microscopic underlying 
theory, it may happen that we have, in an effective theory, properly defined integrals even without it, and that 
no such Hermitisation is needed. This is the case of non-Hermitian Random Matrix Theory with uncoupled 

bosons [laiigin]. 

5.2 Equivalence for real chemical potential 

As already hinted in sect. [2] and explained in sect. 15. II the resolvent method may require the Hermitisation of 
the boson-boson part of Dirac operator : 

/ n"" 2?ei [A + -ri- r A 



~ 1* "/ 



Vet[D, + z,]Vet[D, + ~z,]* , ' 

, nrVet[D. + ,.]Vet[D. + ~z.r ^ + "^^^ , 

lim / m^Vet[D. + z.]VetlD. + ,.r ^ ^ 
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The first line is the proper definition of the resolvent, the second, that is obtained from the first by a trivial 
step, is the Hermitised version, the last line is the regularised integral. Hermitisation is necessary both for the 
computation below and for the underlying microscopic theory, but it's not always necessary for making the 
computation in RMT [iTlITT]. 

In the second line the additional valence quarks are in conjugate pairs {nh+nl\ny+n* + Nf) (the replacement 
Zi — > Zi in the numerator is irrelevant). It will be sufficient to consider the conjugate fermionic quarks as 
independent quarks with mass — m* and chemical potential — concerning the bosonic conjugate quark, as 
already seen in sect. 15.11 there's a deep reason that forces us to treat them simultaneously. 

We mention in passing that a purely bosonic theory with bosons in conjugated pairs (nh + nj, 0) occurs when 
applying the replica trick or Toda lattice equation [9l [161 EH Hi] • 

In the following we will prove the equivalence of the partition functions of exPT and chRMT at non zero 
chemical potential with Uf fermionic quarks with given masses and baryonic potential {mj, /iy), and couples 
of conjugated bosonic quarks, with parameters (mb,Hb) and {—ml,—fj,b)- The partially quenched theory with 
Nf physical quarks, n.^ couples of conjugated fermionic quarks and Uy couples of bosonic quarks will result as 
a special case. The equivalence of the spectra follows applying the resolvent method in both theories. 

We will again use the the two-matrix model as was introduced by Osborn [7] , but the same calculation can 
be done also for the Stephanov model. The equivalent of the QCD Dirac operator for a quark with mass m/ 
and chemical potential /i/ is given by: 



(36) 



where A and B are complex iV_|. x 7V_ random matrices. 

The partition function equivalent to partially quenched QCD is: 



dA dB w{A) w{B) 



XT/ Vet 


ruflN^ iA + jjLfB 
iA^ + lifB^ uiflN^ 






rublN^ iA + HbB 
iA^ + fibB^ niblN- 


Vet 


-mllN+ iA-HbB 
iA''' — ^ibB^ —mllN_ 



(37) 



with a Gaussian weight function 



w 



(X) = Exp [-aNTr [x'' X]] 



(38) 



In order to write the inverse determinants as bosonic Gaussian integrals we perform the same anti-Hermitisa- 
tion as in eq. (j34p . obtaining an anti-Hermitian matrix apart from an e times the identity. The regularised 
denominator is: 



(-) 



= Vet 



Vet 



Vet 



mblN+ iA + fibB 
iA^ + ^.bB^ mf,ljv_ 









£lAr_ 



/e + iA 



V ~^iB^ 



2 

-HB 



rribl 



iA — fibB el 



7V+ 



m+m* 



2 





I 


At - ^fcSt 




iA + fibB\ 











eljv_ 








m+m* 




2 




-iA 




m — m * , 
^ 2 ' 





(39) 



We introduce two sets of iV+ and complex supervectors in (2nb\ni) to write the ratio of determinants 
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as a Gaussian integral 

Zpq = j dA dB d{tp,ip*,(p,(l)*)Exp[-aN {Tr[A^A] +Tr[B^B])] 

X Exp 



(40) 



E 



X Exp 



-E 



'2&-l,/3 



-fibBa/B' 
2 



2 



\ 


1 ^/'2fc-l,a'\ 






026-l,;3' 






V'2&,a' 




J 







^/3a' 2 """/Sq' " 2 

As before we have reduced all the dependence on random matrices in terms of Gaussian functions: 



Exp 



-CtN AaJjA*^^p - iAa,f3 '4'*f,a'Pf,l3 + ^ {'<p2b-l,ahb-l,l3 " '>p2b,ahb,l3) 



"*^",/3 E ^},P'^f^°' + E ('^26-l,/3V'26-l,a " '/'2b,/3'02fc,a) 



(41) 



X Exp 



-aN Ba,fiBl p - Bajj I ^ tJ-fPj^cfkfJJ + E (^2fc-l,a'/'2f,,/3 " '02fc.a'/'26-l,/3) 



and we can perform the Gaussian integration completing the squares: 



Exp 



Exp 



^ 2 



crN 



Q,/3 



a, 13 



(42) 



ctN 



:Str 



a /3 



N N 



Q /3 



Here we have introduced the (2nf,|n^) x (2nb|nj) supermatrices 







and the mass matrix: 







(43) 



(44) 



m/} x Uf^ 

anticipating their TV-dependence below. We can rewrite the partition function ()40p : 



d(V',V^*0, 0*)Exp 



-Str 



N ^ 



(45) 



X Exp 



-Str 



/3 
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From this point on the procedure is the same as before: we integrate exphcitly the sets of supervectors 0, (fi* 
and use the superbosonisation theorem 



Zpq ~ 



d [tp, ip*) Exp 



-Str 



(46) 



X Sdet 



Gl(2nb\nf) 



AT, 



-Str^M-U 



Sdet 
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This result is vahd for finite N . Once performing the N ~f oo weak non-Hcrmiticity hmit f44j, keeping M 
fixed as well as , the following result is obtained 



lim Z. 



pq 



diiH{U)Sdet [C/l^Exp 



Str 



M (c/ + • u^^ ■ r^) + r^r^ur^T^u~^ 



lGl{2nt\nf) 

Rotating all the matrices under the superunitary transformation X ^ ■ X ■ T with 

^/^/2 1/V2\ 
r= I -z/V2 1/V2J 





(47) 



(48) 



we obtain the new result for Hermitised exPT with real 7^ 0, generalising previous results in the literature 
OUT]: 



lim Z„ 



d^iH{U)Sdet Exp 



Gl[2ni\ns) 



Str 



|m ({/ + /• •/) + B^fUB^^^U-^ 



(49) 



where 



M = T'' -M-T 




(50) 



It is expressed in term of S, i^^ and chemical potential using eq. (|3ip . For details on the Haar measure dfiniU) 
we refer to the next section. 

The expression above is a generalisation of exPT with one pair of bosons [24l [17], and of exPT with Nf 
fermions exPT [35]. The signature of the boson-boson block in the metric / differs from [Mj by a factor of 
i because we have used an anti-Hermitisation here, see subsection 15.11 instead of a Hermitisation there. The 
main point is that the two mass terms have different signs if considering bosonic or fermionic quarks, see eqs. 
(102) and (131) of [24], respectively 

The particular cases of eq. (|49p already known in literature [24j \17\ [35] have been derived starting from the 
symmetries of the microscopic theory under vector and axial transformations too, and we suppose that the same 
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arguments can be applied to the general case. The existence of two different matrices Bjf is due to the fact 
that the covariant derivative has a different behaviour on bosonic and fermionic quarks in e^PT as explained 

mm- 



6 Superbosonisation theorem 

The aim of this section is to find a way to express integrals of functions of external products of supervectors in 
terms of integrals over a smaller space. In formulas we are going to study 

J di^lj,r) f(^A®^l^ . (51) 

where i/ife, "0^ are complex vectors oini, bosons and Uf fermions, d{^p, is a shorthand notation for dtpk dip^ 
and / is a function defined over (n;,|n/) supermatrices. 

The main idea is to embed ipk ^ '4'k some manifold where we can define a (5-function and where 
computations are feasible. At a merely symbolic level we have: 

J di^,r) f (^^k®^l^ ^ J^dafia) J di^.r) 6M(^a-Y,i'k®^l^ ■ (52) 

In order to give a meaning to the equation above we have to specify which is the manifold A4 we use and which 
is the measure we use on it (the (5- function depends on it). The use of a J-function requires / to be continuous 
with respect to the metric chosen. 

A good choice is to take as M in eq. ((5^ the manifold of super-Hermitian matrices. One could ask whether 
or not a (5-fmiction may be used dealing with matrices like external product of supervectors [llj : in appendix [Xl 
we show that such a S may be properly defined. The (5-function may be written using an extension of the usual 
Fourier-transform representation of the (5- function on Hcrmitian matrices (see appendix 

A crucial ingredient in this proof will be the possibility to flip the order of integration. Changing the order 
between commuting and anti-commuting variables is not a problem, but problems can arise when considering 
the case of two commuting variables. In order to see where these problems come from, and how to avoid them, 
let's consider a simple example, the one of a single commuting number vector of length 1: 



d^zf{z-z*) = j d^z J dx f{x) 6{x — zz*) 

dxf{x) / d^z5{x-zz*) (53) 
dxf{x) I d^z I dye^^(=^-^^') . 



The first change of variable is always allowed when considering converging integrals. At this point we need 
to flip the order of integration of y and z^z*, and this is an illicit step since the integrals are not converging. 
This problem may be avoided considering a real quantity 77 > 

d^zfizz*) = / d2z/(zz*)e(''-'')^"* 



dxf{x) 



d'z I ^ 






Jm^TT 






' ^ ^^y- f 


ze^ 


-(ri+iy)zz'' 










-1)7: 




t 27r y- 


- irj 

/>OC 




"""Treix) = 


Jo 


dxf{x) . 



(54) 
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and all the steps are mathematically rigorous any time we can apply a (5-distribution as a functional to the 
function f{x)e^^. The symbol 9{x) indicates the step function. 

This easy example is conceptually not too different from the proof of the following theorem. 

Theorem 1 (Superbosonisation) Let f be a function defined on the (n{,|7Z/) supermatrices, then the following 
identity holds 

[ TT d^kd^k fly^A^ i^i] « / d^iHiU) Sdet [Uf f{U) (55) 

J k \ k J JGl(n,\nf) 

whenever the integral on the l.h.s. is well defined. The first integration is performed over N complex supervectors 
in (nfcjny) and dfiniU) denotes the Haar measure over Gl{nb\nf). 

The idea of expressing integrals of a function of an external product ipk ^fpk in terms of an integral of the 
same function over a simpler space has been widely used in physics: it was derived for the first time in [48] for 
external products of Grassmannian vectors appearing in exponential function (the term "bosonisation" comes 
from this kind of application) , subsequently in for commuting number vectors and in [IHl ISHl Ell IMl E3 
for supervectors. This (super)bosonisation may be seen as application of the Riesz Representation Theorerro- 
However, up to our knowledge, a graded version of the last theorem is lacking. 

By coincidence the superbosonisation theorem has been independently developed at the very same timtd 
of an analogous theorem on superbosonisation by Littelmann, Sommers and Zirnbauer [101 111] . Anyway the 
proof they give is different from ours, theirs is an algebraic proof, ours uses analysis instruments. The main 
difference is that we base our proof on the existence of a i5-function fulfilling eq. ([5^. and that collaboration 
develops a powerful apparatus in order not to use such an equation fTUI . For this reason we have chosen to show 
in appendix E] all the details concerning the mathematical rigorousness of our definition. 

Proof: As already hinted above we write the l.h.s. of eq. introducing an additional integration over 

supcr-Hcrmitian matrices: 



fe 



H=Hf 




OC / di//(i/)e''«*'-[^l / d(^,V^*)e-"5:.V'I--^. / dFe'^'i"-^^^'^^'") 
Jh=h^ J Jf=f^ 

= f dH f{H)e''^'''^"^ f dFe'^" Sdet[ri + iFy^ (56) 

JH=H'< J F=F'< 

where d(V', V'*) = 11^ dtpkdipl and dtpkdtpl = JJ'^'' d^tpk,b H/^ dipkjdipl j. We denote by = -0*'^, and 

is the external product and ip^ ■ ip the scalar product in (nf,|ny). The measure dH is the flat measure in the 

independent entries of the super-Hermitian matrix. 

From now on let's focus on the second integral. This quantity reminds of a quantity already computed in 
|53] . unfortunately we cannot use this result because it does not take care about the boundary terms arising 
from the diagonalisation of super-Hermitian matrices (Efetov-Wegner terms). In the following we will never 
change the Z-gradings of the integration manifold, and hence no boundary terms will arise [55) . 

Writing F and H in terms of blocks, dF = dFi dF2 d^ d$t 

Fi $t\ et 

we obtain: 



I = dFe'^" Sdet[F -iif^-^ (58) 

J F=F'< 

= J dFi dF2 d^ d$t / 1 1 Exp [iTr [F^Hi + $^9] - iTr [H2F2 + $6^] ] . 



°This theorem ensures that a bounded Hnear functional of functions over a locally compact Hausdorff space may be computed 
like an integration of that function over that space using a proper measure [54] . 

''This work was already presented by one of the authors in "QCD in extreme conditions" , Frascati 6-8 August 2007. 
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As a first step we note that if we diagonalise tlie matrix Fi and make an analytic continuation in tlie 
eigenvalues we see that the only poles in the expression above are the ones where at least one eigenvalue of 
Fi is irj. Considering the integration on the eigenvalues as a complex contour integral and closing the integral 
in the upper semicircle (or lower one, depending on the signature of the matrix Hi H) we have non vanishing 
contributions only if every path of the eigenvalues winds the pole in ir]. We can consider this integral as a 
contour integral around irj whenever 9{Hi) > 0, otherwise it's zero. Since there are no poles apart from the 
ones in ii] we can fix the contour integral as we prefer. We choose the modulus of the eigenvalues as equal to 
1, hence we have that the matrix Fi — irj may be analytically continued obtaining a unitary matrix. 



N 



X Exp [i [FiHi + $te] - i [H2F2 + $et]] . (59) 

The notation fjj^^^^-^ dFi stands for an integral over the manifold U{ni,) considering as a measure the analytic 
continuation of the flat measure of Hermitian matrices. It is related to the Haar measure c?/i[/ integration by: 

dFi^ [ dfiu{Fi)Vet[Fi]"'' . (60) 

The relation between these two integrals is the same as between the two integrals below describing a circui- 
tation around zero, with z — 



ie. 



27r /. 2 



d6l = * dz- ^ f dz. (61) 

J\z\ = l ^ J\A=^ 

Eq. ((50)) can be derived diagonalising the matrix Fi = U ■ f ■ W where fj = e*^^ : 
/ dFi = TT i"rff,A2({f}) / dU 



2?et[f] A({f})A({ri})l?et[f]--i 

i 

d9,\A{{e'''})\^Vet[f\'"' J dU 



dU 



We consider now that 



d^iu{Fl)Vet[Flp . (62) 



-iN -iTr[{A-i'n)H2] 



f dAVet[Af e-'^"^^^"'^ = / dA Vet [A - ir,f e- 
Ja=a'i Ja=a'i 

= I dAI?ei[A-z77-<i>tFfi$]'^e-'^'-M-"'-*^i"*')^=] . (63) 
Ja=ai 

The first equality comes from the analyticity of the integrand in the diagonal entries of the matrix and 
the second comes from the analogue of the contour invariance in superanalysis [56^ . applied to the real and 
imaginary parts of the Hermitian matrix entries. Applying this equivalence we obtain 



I = 9{Hi)e 



■Juin,,) J 



N 



U(n,) J \Vet[Fi] 
Exp [iTr [FiHi + $^6] - iTr [H2F2 + $6^] - iTr [i/a^i^f^*^]] ■ (64) 



°Thc matrices Hi and H2 are complex number Hermitian matrices. 
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The subsequent step is to perform the $, $^ integration, using 

Tr [$^e] - Tr [$6^ + i/a^i^f^*^] = -Tr [9$^ + $6^ + Hs^-Ff^^'^] 

= -Tr [OFiFf + H2^Q^H2^ + H2^F^^<i>^ 

= -Tr[{eFi+H2<i>)- {F^^<i>'' + e^H^^) -QFiQ^'H^^]. (65) 

We can transform the fermionic variables of integration $ ^ $ = + QFi and $t ^ $t ^ F^^^^ + 
Q^H:^^ and thus perform the Gaussian integration, 



d$d$tExp [-iTr [i/a^-Ff^*^] + iTr [QFiQ^ H:^^]] = 



Vet [ffz]"' ^iTr[efietg-M 
Pet [Fi 



1"/ 



t -iTr[*l>+l 



where the last integration is just a constant. Getting back to eq. we have: 

T oc e{H,)Vet[H2r e-^'^'^i^U dF, [ dF2 ^^^^^^ g.Tr[F,g,]-»T.[g.f.]+.Tr[eF,etg-^] 



(66) 



(67) 



= 6{Hi)Vet[H2T'' e 



rib -riStr[H] 



dF2Vet [F2Ye 



N -iTrlH-iFn' 



dFi- 



Uiub) ^Vet[F,f+''' 



^'i\Fi{Hi-e'<H-^e)]_ 



Let us consider again the whole eq. (|56p . We can now perform the integration in Fi and using the 
integrals in appendix [Bl 

/d(V',V'*)/( VV'fc^V'i ) oc / dHi9{Hi) f dH2Vet[H2T'' ( dQdQ^ f{H) 
J \ J, / Jhi=h\ Jh2=hI J 

X / dF2Vet[F2f e-'^'-^''-^-^ i dF, L^e4^i(«i-etff.-^e)] 

Jf2=f^ Jum Vet[F,f+^' 

oc / dHxe{Hi) I diJa X>et [iJa]"" [dede^f{H) 

X / dF2 Vet [F2f e-'^'-I^^^^lpei [H, - H^^Q] ^+"^-"' 



dHi6{Hi) 6 dH2 I dSde 



^ Vet [Hi - e'<H^^e] 



— 1 N+nj-rib 



Hi=Hl 



dHi9{Hi) i dH2 / dQdQ'^Sdet 



U{ni) 



Vet [i72]'^+"^'""'' 

'Hi et 



-f{H) 



e H2 



N+nf-Tib 



f{H) . (68) 



The last step missing is to find a relation between the measure used in the equation above and the Haar 
measure. This relation is: 



Hi=h\ 



dHiOiHi) f dH2 I dQdO^Sdet 



Hi et 
e H2 



- Uf-Ub 



dflH 



Gl{nb\nf) 



Hi et 
e H2 



(69) 



It is proved in the following. The manifolds we are integrating are the same, in fact the Hi = Hi > is equivalent 
to the boson-boson [40j base manifold of Gl{nb\nf) that is Gl{nt,) /U{ni,) [421I29]. The analytic continuation of 
the eigenvalues of H2 = H\ making a circuitation around zero is equivalent to the fermion-fermion part U{nf) 

HZ]- 

We have now to consider the integration measures in eq. (|69p . The point is that the measure on the l.h.s. 
is the Haar measure of super-Hermitian matrices (a group closed under addition) computed on the analytically 
continued manifold. It is induced by the flat metric: 



Str[dU-dU] , 



(70) 
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but the Haar measure on the r.h.s. is the analytical continuation of the Haar measure on super-Unitary matrices 
[m (a group closed under multiplication) , induced by the metric |^ : 



Str [dU ■ dU-^] = Str [U-^dU ■ U-^dU] 



(71) 



In order to find a relation between these two matrices one has to compute the Berezinean of the transfor- 
mation [SI [12]: 

'{U-^dU)i 



05 = SdeU 



dU„ 



(72) 



Since this is not a change of variables (the differentials dU are the same) and hence no Z-gradings are 
changed, no boundary (Efetov-Wegner) terms arise in the superintegration. We write this Berezinean matrix 
in the block structure (Latin indices stay for bosons and Greek for fermions): 



U-^dU 
dU 



/bb/bb bb/ff bb/fb bb/bf\ 

ff/bb ff/ff ff/fb ff/bf 

fb/bb fb/ff fb/fb fb/bf 

\bf/bb bf/ff bf/fb bf/bfj 



/Kb' ® In. 

Uj ® 1„, 



Ua0 ® 1. 







V uj 1„, 

Computing the superdeterminant of the matrix above we obtain 

'{U-'^dU), 



Uj ® In, 



(73) 



As a consequence we have: 
d^iH 



<B = SdeU 



Hi et 
e H2 



dU„ 



Sdet [U] 



rif—nt, 



dHi dH2 de Sdet 



e H2 



(74) 



(75) 



that is eq. l[^ . Together with eq. ([^ we obtain the proof of our theorem eq. l[S5|) . 



7 Conclusions 

In the present paper we prove that xPT in the e-regime and chRMT have the very same spectral properties 
in the microscopic (weak non-Hermiticity) limit. The equivalence holds for zero, imaginary and real chemical 
potentials. This result is achieved by proving the equivalence (up to an irrelevant overall constant) of the 
partition function of the partially quenched chRMT, with bosonic and nj fermionic quarks, to an integral 
defined over the maximal Riemannian subspace of the linear group of {nb\nf) supermatrices for any finite N. 
Performing the N oo limit of this integral one obtains the partially quenched e^PT. The one obtained for real 
chemical potential was derived here for the first time and agrees with the previously known purely fermionic and 
purely bosonic cases. An essential ingredient of this computation is a superbosonisation theorem connecting 
integrals of external products of supervectors to integrals over Gl{ni,\nf). An explicit parametrisation of integrals 
over Gl{nb\nf) is provided as well. 

The equivalence we have established by matching the generating functionals of density correlations holds 
also in cases where the densities are not yet known explicitly. Only for at most two different chemical potentials 
all densities have been derived so far. Our results also shed some light on universality in chRMT with chemical 
potential by matching the Gaussian one- and two-matrix model of Stephanov and Osborn, respectively. This 
gives hope to establish a further reaching universality beyond Gaussian RMT. 
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A Defining (5-functions on super-Hermitian matrices 



The question is whether or not one can define a ^-distribution such that 

f(^i'k® - J dM f{M)5 (^M - V'fe ® V'i^ • (76) 

In order to clarify the idea we take the same simple but nontrivial example as in |1 Ij . that is k = 1, 
{nh\nf) — (0|2). The external product is a 2 x 2 matrix, it's terms are only nilpotent commuting numbers 
and -0* for i, j = 1,2. We can consider this matrix as belonging to a superanalytic continuation of a manifold 
of real dimension 4 (at least). We could take as a such manifold C/(2), Gl{2)/U{2), Hermitian matrices or real 
matrices. Despite in the rest of the work we have used Hermitian matrices, for this simple example we will use 
the one with the simplest notation, real matrices: 



a : 



Once considering the external product of vectors like an element of a commuting number manifold, the well 
definiteness of the ^-distribution is inherited from the one of the base manifold by superanalyticity in the real 
parameters of the manifold. In formulas: 

/(a) EE /(aii,ai2,a2i,a22) = y dm (5(7ti - a) /(mn, mi2, m2i, 77122) , (78) 

where the integration manifold is the base manifold (real matrices) J dm = J^^Yli j=i dinriij, and the 5 distri- 
bution is the one of the base manifold too, 5{m — a) = Jl? j=i ^{''^ij ~ Qy)- 

The last thing to be verified is that the superanalytic continuation of the ^-distribution is properly defined; 
this was already shown in [55]. Given a quantity x = xb + xs, where the first part is the body of the number, 
and the second is the nilpotent part. Given a function with a sufficient number of derivatives in xb such that 
f{xB+xs)^ I]r?=o /^"''(^s)^ is properly defined, we have: 

dy6{y-XB-xs)f{y) ^ / dy 5^-\y ^ x bY-^ f [y] 
-00 J-00 n. 

E^^/ dyS{y-XB)f^''\xB) = f{xB+xs). (79) 



n=0 



Completing the discussion on the example above we can write explicitly d{m — a): 

2 

S{m - a) = 5{m) - ^ 'ipiip*dm^.S{Tn) + V'iV'i'02V'2(c^miif^m22 + dyyi^^dyn^^)5{m) (80) 
where m is a real number 2x2 matrix. 



The generalisation for what is done above to eq. (|76p is straightforward: one considers the external 
product as an element belonging to the superanalytic continuation of a manifold. The manifold we choose is 
the one of super-Hermitian matrices: the boson-boson and the fermion-fermion block are Hermitian matrices, 
and the 5 is defined as the superanalytic continuation of the product of (5's over the real and imaginary parts 
of the independent elements of the matrices; boson-fermion and fermion-boson blocks are made of independent 
Grassmann integration variables and 5 may be represented as: 

S{d-d)<x{e~0)<x f d^e'i'-''-^\ (81) 
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As an exponential representation of S is allowed both for commuting and anticommuting variables, such a 
representation of (5-function on super-Hermitian matrices is allowed too: 



/(«) 



Af=Aft 



dM f{M)S{M ~ a) 



(82) 



dM f{M) 



^iStr[F{M-a)] 



F=F1 



where dM and dF stand for the flat measure over the independent entries of super-Hermitian matrices: 



dM = JIJI dM^^a dMa,i Jl dM,,j J| dRe ^ dim M,j Y[ dMa,a Y[ ^'^^ ^^.p dim 



/3 



(83) 



i—l Q— 1 



a, (3 



B The integrals in eq. (1671) 

B.l Boson-boson block 

Let's start from 



dFi 



1 



Vet [F,^^+''f 



Exp 



iTr 



FiH 



diiu[F^)Vet [^1]-^-"^+"" Exp 



(7(rii, 



iTr 



FiH 



Vet 



H 



N+Uf—nj, 



dnu{Fi)Vet 



FiH 



Exp 



iTr 



FiH 



(84) 



where for simplicity of notation we use H = Hi - ^6^- The integral in the last line may be performed 
using the character expansion. This is a particular case of a more general integral computed in [57]; the result 
is: 



(0) 



(85) 



where the sum is over the irreducible representations of Gl{ni,). The quantity Xr(0) is zero for all the represen- 
tation apart from the trivial one. The result is just a constant 



dFi 



Vet [Fif+''f 



H 



N+rif-ni, 



(86) 



This result, together with the argument of analytic continuation of Hermitian matrices, gives an alternative 
way for computing the Ingham-Sicgcl integral as was done in [27j . 

B.2 Fermion-fermion block 

The fermion-fermion block integral has already been computed by Guhr [53 . The result is: 



/ dF2Vet[F2f Exp[-iTr[F2H2]] oc TT (/i^- 
Jfo^fI 



(87) 



Here we don't have any problem with possible Efetov-Wegner terms arising from these diagonalisation: both 
F2 and H2 are complex Hermitian matrices and no Z-grading will be changed during diagonalisation |55| . 



^In order to compute the exact value one has just to substitute the relation a 
are the labels of the representation. 



(") 



detw 



, where r = (m , . . . , n^) 
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If we integrate this functional with a given function analytic in the matrix entries, using the Cauchy integral 
formula we obtain: 



dH2 
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